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ABSTRACT: With the large-scale integration of renewable en- 
ergy and grid-following power electronic devices, system volt- 
age support strength decreases, thus threatening the system's safe 
and stable operation. In the homogeneous scenario where power 
electronic devices integrated into the system have the same dy- 
namics, a theoretically rigorous and highly stable quantitative 
method for quantifying system voltage support strength can be 
formed based on the generalized short-circuit ratio(gSCR) and 
the device or station critical short-circuit ratio(SCRo) In the het- 
erogeneous scenario where power electronic devices integrated 
into the system have weakly different dynamics, system voltage 
support strength can be quantified by the first-order approxima- 
tion of gSCR and SCRo, based on some special characteristics of 
devices and power grid. However, there is a lack of unified der- 
ivation principles and calculation methods. To this end, this pa- 
per focuses on quantifying voltage support strength under small- 
signal stability and discovers that the multi-infeed system can be 
approximately decoupled into multiple low-dimensional sys- 
tems. On this basis, this paper proposes the concept of eigen- 
subsystem and its calculation method and interprets their physi- 
cal significance. Based on the concept of eigen-subsystems, the 
derivation principle and calculation method of gSCR and SCRo 
are given. Additionally, the specific calculation methods for in- 
dices are provided when considering grid-following converters 
under non-rated operating conditions, reverse active power out- 
put and grid-forming devices. Finally, the effectiveness of the 
principles and methods is verified in several cases. 
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摘要 : 随 着 跟 网 型 新 能 源 /电力 电子 装备 的 大 规模 接 入 ， 系 
统 电压 支撑 强度 降低 ， 新 型 电力 系统 的 安全 稳定 风险 增加 。 
在 同 构 新 能 源 /电力 电子 装备 接 入 场景 下 ， 利 用 电网 广义 短 
路 比 和 装备 / 场 站 临界 短路 比 可 以 形成 与 稳定 性 强 相 关 的 系 
统 电压 支撑 强度 量化 方法 ; 在 弱 异 构 场 景 下 , 基于 装备 和 电 
网 动态 的 特殊 性 质 ， 利 用 广义 短路 比 和 装备 / 场 站 临界 短路 
比 的 一 阶 近似 可 以 实现 强度 量化 , 但 缺乏 统一 的 导出 原理 及 
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计算 方法 。 为 此 , 本 文 聚焦 小 扰动 维度 下 电压 支撑 强度 量化 
问题 , 首先 利用 多 馈 入 系统 可 近似 解 耦 为 多 个 低 维 子 系统 的 
规律 , 提出 特征 子 系统 的 概念 和 计算 方法 , 并 诠释 了 其 物理 
意义 ; 其 次 , 基于 特征 子 系统 , 提出 了 广义 短路 比 及 装备 / 场 
站 临界 短路 比 的 导出 原理 及 计算 方法 。 此 外 , 针对 电力 电子 
装备 在 非 额 定 运行 点 、 部 分 装备 有 功 功 率 反 向 以 及 包含 构 网 
型 装备 的 场景 下 ,给 出 了 广义 短路 比 的 具体 计算 公式 。 最 后 ， 
算 例 验 证 了 所 提 原 理 和 方法 的 有 效 性 。 


关键 词 : 广义 短路 比 ， 装 备 / 场 站 临界 短路 比 ， 特 征 子 系统 
0 引言 


随 着 新 能 源 的 快速 发 展 ， 电 力 系 统 逐 步 演化 成 
含有 高 比例 新 能 源 及 电力 电子 装备 (后 文 简称 为 
“装备 ”) 的 新 型 电力 系统 中 中。 电压 支撑 强度 表征 
了 扰动 下 系统 电压 抵抗 偏 移 以 及 抵抗 失 稳 的 能 力 ， 
强度 高 表明 系统 抵御 被 外 力 破坏 的 能 力 强 ， 受 扰 后 
电压 响应 性 能 好 。 然 而 ， 具 有 强 支 撑 作 用 的 同步 机 
占 比 减少 ， 新 能 源 集 中 并 网 容量 大 、 与 主 网 电气 距 
离 远 且 常 采用 跟 网 型 控制 策略 而 对 电网 支撑 强度 
弱 ， 使 得 新 型 电力 系统 电压 支撑 强度 变 低 ， 受 扰 后 
安全 稳定 风险 增加 B77， 吸 需 科 学 的 强度 量化 和 计 
算 方法 (考虑 到 低 惯 量 维度 的 强度 问题 可 以 分 开 考 
谍 ， 后 文 为 表述 方便 ,将 系统 电压 支撑 强度 简写 为 
“系统 强度 ”)。 

大 规模 跟 网 型 装备 (如 不 特别 说 明 ， 后 文 简称 
“装备 ”) 接 入 后 , 系统 强度 的 分 析 计 算 存 在 以 下 难 
点 : 装备 及 电网 母线 数量 多 ， 且 装备 与 网 络 动态 交 
互 复杂 ， 高 维 动态 特征 难以 简化 ， 装 备 种 类 多 ， 如 
空 制 策 略 和 参数 等 各 异 的 跟 网 型 变 流 器 或 构 网 型 
变 流 器 等 BTINl， 进 一 步 增加 了 系统 的 复杂 度 ; 新 能 
源 运行 工 况 多 ， 也 给 强度 量化 带 来 巨大 挑战 。 

为 分 析 系 统 强 度 且 方便 工程 应 用 ， 工 业界 常用 
电网 强度 近似 表征 系统 强度 ， 并 引入 单 馈 入 和 多 人 馈 
入 短路 比 作 为 量化 指标 。 其 中 ， 电 网 强度 表征 了 在 
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不 计 及 跟 网 型 装备 动态 作用 下 母线 电压 在 扰动 作 
下 抵抗 偏离 正常 运行 点 的 能 力 。 例如 : CIGRE 多 
馈 入 短路 比 (multi-infeed short-circuit ratio, 
MSCR)H31, 位置 相 关 短 路 比 (site-dependent SCR, 
SDSCR)!!41, 4 Yq wh Fae LL (multi-renewable energy 
station SCR, MRSCRJ)U5、 等 效 短路 比 (equivalent cir- 
cuit-based SCR，ESCRJ)023 和 带 交 互 因子 的 短路 比 
(the SCR with interaction factors, SCRIFP)D] 等 指标 。 
上 述 指标 可 以 反映 装备 接 入 前 后 母线 电压 稳 态 1 
移 量 ,但 由 于 未 考虑 装备 间 及 其 与 电网 的 宽频 耦合 
特性 ， 故 与 系统 稳定 性 关联 弱 , 用 于 反映 稳定 特性 
I oR ZS LER SCTE 

文献 [7]、[16]-[18] 从 装备 和 电网 动态 交互 视角 
E 塑 了 融合 装备 和 电网 动态 的 系统 强度 内 涵 ， 并 提 
基于 电压 对 电流 (或 电压 对 功率 ) 灵敏 度 方程 的 
电网 广义 短路 比 (generalized short-circuit ratio, gSCR) 
指标 和 装备 / 场 站 临界 值 指标 , 从 理论 上 实现 了 系统 
强度 与 小 扰动 稳定 性 的 关联 。 在 同 构 系 统 〈 即 接 入 
装备 的 动态 特性 相似 ， 有 具体 定义 见 文 献 [7]) 中， 多 
电力 电子 装备 馈 入 系统 (简称 “多 馈 入 系统 ”) 可 解 
耦 为 装备 动态 整体 保持 的 多 个 等 效 单 馈 入 系统 ; 基 
等 效 单 馈 入 系统 ， 类 比 传统 短路 比 的 用 途 ， 定 义 
电网 广义 短路 比 和 装备 / 场 站 临界 短路 比 , 并 论证 
了 装备 / 场 站 临界 短路 比 等 于 广义 短路 比 临界 值 的 
寺 论 ， 从 而 实现 了 考虑 装备 动态 特性 的 系统 强度 量 
; 在 弱 异 构 场 景 下 ， 如 外 特性 存在 差异 的 直 驱 风 
机 上 、 光伏 本 、 构 网 型 装备 0000 和 储 能 系统 接 入 09] 
或 者 装备 在 非 额定 运行 点 [等 场景 , 异 构 装备 动态 
可 看 作 是 同 构 系统 中 装备 或 网 络 动 态 的 扰动 ， 利 用 
广义 短路 比 和 装备 / 场 站 临界 短路 比 的 一 阶 近似 即 
可 实现 强度 量化 。 上 述 广 义 短路 比 及 其 临界 值 的 计 
算 关 键 为 求解 能 够 近似 原 系统 关键 模 态 稳定 性 的 
解 耦 单 馈 入 系统 。 然 而 ， 在 这 些 非 理想 场景 中 ， 寻 
找 等 效 解 耦 单 馈 入 系统 的 理论 和 方法 缺乏 普 适 性 ， 
不 利于 广义 短路 比 在 复杂 场景 下 的 应 用 。 

为 此 , 本 文 则 在 统一 多 种 场景 下 广义 短路 比 的 
导出 原理 及 计算 方法 。 首先 , 引入 特征 子 系统 的 概 
念 并 提出 其 求解 方法 ; 其 次 , 基于 特征 子 系统 提出 
广义 短路 比 与 装备 / 场 站 临界 短路 比 的 通用 导出 原 
理 和 计算 方法 ; 最 后 , 算 例 验证 了 方法 的 有 效 性 。 
1 预备 知识 
1.1 FERRET SST CREA 

2nx2n FERE M 的 特征 子 空间 及 广义 特征 矩阵 

定义 如 下 : 

定义 120〈 特 征 子 空间 和 广义 特征 矩阵 ): 记 
R(X) ={Xu:u € C”} FRY) = [Yuu eC") 分别 为 


=r 


m. cr im 
DE 
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ERE M 的 右 和 左 特征 子 空间 , Sah X YAA 
量 张 成 , 归 一 化 后 满足 YX=1; 与 特征 子 空间 对 应 ， 
满足 式 () 的 工 和 矩阵 称 之 为 矩阵 M 的 “广义 特征 矩 
阵 ”， 其中， 广义 特征 矩阵 的 特征 值 也 是 矩阵 M 的 
特征 值 。 


MX = XL, Y*M = LY" 
由 上 式 (1) 可 知 , L 可 由 下 式 求解 得 到 : 
L=Y"MX 

广义 特征 和 矩阵 及 特征 子 空间 分 别 是 特征 值 及 
句 量 的 推广 ， 满足: 若 (v, DE REEE L 
的 特征 对 ( 式 中 * X L ARIE EE), MWA, AN 
JERE M 的 特征 对 。 其 中 ， 特 征 值 是 维 数 为 1 的 特 
殊 广义 特征 和 矩阵， 特征 向 量 张 成 的 空间 是 特征 子 
空间 的 特例 ， 对 应 L 为 对 角 和 矩阵 。 为 方便 理解 及 
后 文 推导 ， 下 文 举例 说 明 。 

假设 矩阵 和 MM 有 一 组 特征 子 空间 R(X,) 、R(Y) 
(二 1,...,n)， 对 应 有 个 2 维 广义 特征 矩阵 五 ， 表 达 
式 如 式 (3) 所 示 ， s X=(Xa, X2) Y=(Yn, Yn)» 
Xas Xos Yas YoER2""!. 根据 特征 子 空间 的 定义 
可 知 ， 利 用 和 矩阵 CQ ... XDM: ... YAK M 转换 
为 块 对 角 抢 阵 ， 且 对 角 块 为 广义 特征 矩阵 Lo 表示 
为 : 


(1) 


(2) 


特征 


(Y,...Y,)" M(X,...X,) = diag(L) 


式 中 : Y” MX, =L 
由 于 (Xi ... X,) AIM... YHET E E AE 
且 列 向 量 分 别 线性 无 关 , 故 式 (3) 的 矩阵 变换 不 改变 
JERE M 的 特征 值 ， 即 矩阵 diag( 五 ) 与 M 的 特征 值 
完全 相同 ， 其 中 Li 的 特征 值 为 M 的 部 分 特征 值 ， 
Li 可 理解 为 M 经 过 特征 分 解 得 到 的 二 个 解 簿 子 矩 
o Ib, Æ MA 个 2 重 特征 值 Li=1,.…….,n) 
且 存在 两 个 线性 无 关 的 特征 向 量 ， 那么 其 张 成 的 空 
间 为 特征 子 空间 ， 此 时 广义 特征 矩阵 满足 Lho 
1.2 ”特征 子 空间 摄 动 
考虑 如 下 受 扰 矩阵 : 
M=M+E 
A, ERE E NIERE M 所 受 扰动 。 
若 扰动 矩阵 E 足够 小 ， 用 的 特征 子 空 
的 特征 子 空间 存在 近似 关系 ， 如 下 引 理 。 
引 理 120 (特征 子 空间 摄 动 ): TEE M 的 右 
和 左 特征 子 空间 左右 分 别 为 只 和 ) MRY), WA 
wr REEF 28 ta] R(X) AM R(F) 以 及 广义 特征 矩阵 
的 近似 满足 : 
R(X) = R(X +o(|EJ)) RL )= R(X +o(]E])) © 
(6) 


(3) 


(4) 


= a] M 


TH 


L=F" MX =Y" MX +ol(|E| ) 
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其 中 , o(.) 表 示 高 阶 无 穷 小 , | 赂 表示 矩阵 任意 范 数 。 

由 引 理 1 可 知 , M 的 特征 子 空间 可 以 近似 让 的 
特征 子 空 间 。 扰动 越 小 , 它们 的 特征 子 空间 越 接近 。 
当 六 的 特征 子 空间 计算 非常 复杂 时 ， 可 构造 合适 的 
M 简化 下 的 广义 特征 矩阵 计算 。 附 录 A 给 出 了 两 
维 矩 阵 的 例子 说 明 上 述 结论 。 
1.3 ”特征 子 系统 

考虑 2n 维 的 动态 系统 4， 记 闭环 传递 函数 矩 
阵 和 系统 矩阵 分 别 为 TAOM MTS), ASAE TT 
程 det(M)=0 决定 了 A 的 稳定 性 。 根 据 定 义 1 并 
类 比 式 (3) 可 知 , 站 有 一 组 特征 子 空 间 和 多 个 低 维 广 
SURE REL, AL 的 特征 函数 也 是 让 的 特征 函 
数 。 因此， 广义 特征 矩阵 上 (ke 1,... RPE R 
数 能 够 决定 系统 在 某 一 关注 模 态 下 的 稳定 性 。 为 
利用 低 维 的 工分 析 用 矩阵 的 性 质 ， 引 入 了 特征 子 
系统 的 概念 ， 具 体 为 : 


4， 其 系统 闭环 传递 函数 道 和 矩阵 称 之 为 该 系统 的 系 
BERE, BN =T(s); 记 必 的 第 个 特征 子 空间 和 广 
REEERE A K(X). RP) ML, RER 
M) RREY L 的 低 维 子 系统 称 之 为 4 的 特 
征 子 系统 。 


L =Y,'MX,,k e{l pins n} (7) 
2 系统 建 模 及 问题 描述 
设备 类 型 多 样 运行 方式 多 变 


首先 ， 考 虑 所 有 线路 阻抗 比 相同 ， 从 装备 端口 
往 交 流 电 网 看 ， 在 全 局 同步 坐标 系 下 的 电网 端口 电 
压 微 增 量 AUww(s) 和 电流 微 增 量 ALw(s) 的 灵敏 度 方程 
可 表示 为 L191: 


-AL =(B@6G(s))AU,, (8) 
Ju (9) 
Bı By +m+r) 
WH: B=| : : 
Bl . an 
人 ee ee 1 | 
s*|(aq+ty +1, -l s/(@ +7) 


其 中 , B 为 包含 所 有 节点 的 等 效 导 纳 矩阵 或 灵敏 度 
矩阵 ; r=R/X 为 线路 阻抗 比 ; A 表示 变量 的 增 量 ，@® 
表示 Kronecker FR; wo 为 交流 电网 的 同步 转速 。 当 
网 络 侧 包含 电力 电子 无 功 补偿 、 构 网 型 装备 或 负荷 
等 元 件 时 ， 修 改 交 流 网 络 端口 电压 对 电流 的 灵敏 度 
矩阵 即 可 考虑 这 些 元 件 带 来 的 影响 ， 后 续 强 度量 化 
方法 仍然 适用 ， 限 于 篇 幅 不 再 详细 展开 。 
其 次 ,装备 端口 的 电压 微 增 量 AUww(s) 和 电流 微 
增 量 ALy(s) 的 灵敏 度 方程 可 表示 为 (9: 
Yee (s) 


Yipe (8) 
EP, Vieri(s)~ Yerom (ONIA ntm 台 装 备 的 导 
纳 传递 函数 矩阵 ，02>ER2”“2…， 元 素 均 为 0。 
最 后 , 结合 式 (8) 和 (9), 利用 如 下 行列 式 可 表示 
图 1 所 示 系 统 的 闭环 特征 方程 ; 
det(T,,, (s)) =det(Yine(8)+Y¥u(s))=0 (10) 
A, Tres) AE AY AIA Fe eb Pe A FE ME, BY 


1 多 类 型 新 能 源 馈 入 电力 系统 
Fig.1 The new power system with multiple types of renewa- 
ble energy devices integrated 


考虑 图 1 所 示 新 型 电力 系统 ， 包 含 ntm RX 
备 及 接 入 的 母线 ，r 个 中 间 无 源 母线 和 1 个 无 穷 大 
母线 ( 若 有 构 网 装备 , 无 穷 大 母线 可 以 没有 ), 其 中 ， 
n 台 装 备 是 新 能 源 或 储 能 的 跟 网 型 变 流 器 ; m 台 是 
构 网 型 装备 ， 如 构 网 型 风机 、 储 能 或 调 相机 等 。 
2.1 系统 动态 建 模 

在 全 局 xy 坐标 系 下 ， 建 立 图 1 所 示 系 统 的 线 
性 化 模型 ， 得 到 系统 的 闭环 特征 方程 。 


ERRER, de(l) KR RITI. 
22 ”基于 广义 短路 比 的 系统 强度 量化 思路 回顾 

式 (10) 所 示 闭 环 特征 方 程 决定 了 系统 稳定 性 ， 
然而 ， 由 于 Tyas) 维 数 较 高 ， 故 寻找 近似 的 低 维系 
统 以 简化 分 析 。 

对 于 同 构 系统 ， 将 多 馈 入 系统 解 耦 为 多 个 等 效 
单 馈 入 系统 。 基 于 失 稳 风险 最 高 的 等 效 单 馈 入 系统 ， 
类 比 传统 短路 比 定 义 ， 可 给 出 原 系 统 的 广义 短路 比 
gSCR 与 装备 / 场 站 临界 短路 比 SCRo, 并 形成 基于 广 
义 短 路 比 的 系统 强度 分 析 方 法 ， 而 且 具 有 电网 动态 
与 装备 动态 分 离 的 特点 。 计 算 公 式 如 下 式 所 示 吕 。 

gSCR > CgSCR, = SCR, (11) 
g a 8SCR-SCR, 

SC 


或 B x100% > 8, % (12) 


4 


ERN: 


aL zl a l 
gSCR = An Sint red Soins) =o"! (s Bimm) Z red S Blasm) ( 1 3) 


式 中 : Saw+m =diag(Sai,..., Saintm), diag) KRX FA FB 
ME, Sei., Sptm) 分 别 为 ntm 台 装 备 的 额定 容量 ; 
Brea NR FA AS Pr PEER AY PH 2 A E BE ER EB As 


中 E HL 


FEA, CgSCRo 为 广义 短路 比 临界 值 ; fo% 是 为 了 保 
证 系统 具有 足够 的 裕 度 而 设置 的 阔 值 ，gSCR 的 表 


X — 
BILK 


BEERE; Zer Bed 5(:) 和 4 O 分 别 表示 求解 


和 矩阵 的 最 大 奇异 值 和 最 小 特征 值 。 


广义 短路 比 是 灵敏 度 矩阵 的 最 大 奇异 值 ， 反 映 
了 多 端口 交流 网 络 中 电压 扰动 对 电流 扰动 的 最 大 


灵敏 度 [7。 


>4 


于 异 构 系 统 ， 利 用 和 矩阵 摄 动 理论 对 gSCR 和 


SCRo 进行 近似 即 可 。 具体 地 , 根据 摄 动 来 源 以 及 装 
备 和 电网 动态 特性 分 为 两 类 : 1) 跟 网 装备 动态 : 跟 
网 型 变 流 器 外 特性 存在 差异 ， 认 为 对 同 构 系 统 装备 


动态 存在 摄 动 ， 


此 时 计算 gSCR 的 边界 条 件 不 变 ， 


而 SCRo 取决 于 装备 动态 特性 的 加 权 平 均 ， 具 体 见 
文献 [7][17][19][20];” 2) 等 效 交 流 网 络 动 态 ， 当 考 


es 


虚构 网 型 变 流 器 或 调 相 机 等 电压 源 特性 占 优 的 接 
入 装备 ， 可 认为 是 对 网 络 修正 ， 此 时 SCRo 的 计生 


方法 不 变 ， 而 修正 gSCR 的 计算 边界 条 件 即 可 计 及 
这 些 构 网 装备 的 影响 ， 有 具体 见 文 献 [9][10][11]。 


2.3 拟 解决 的 问题 


综 上 所 述 ， 


gSCR 指标 及 其 临界 值 的 定义 和 计 


算 ， 其 关键 步骤 是 寻找 等 效 的 低 维 单 馈 入 系统 ， 即 
特征 子 系统 。 然 而 ， 现 有 研究 仅 适 用 于 单一 典型 场 


景 下 的 问题 ， 如 


[9、 构 网 型 装备 


仅 考 虑 外 特性 异 构 的 跟 网 型 变 流 器 
1000 或 储 能 系统 0 等 装备 接 入 的 场 


AI» 其 方 法 及 理论 


论 缺 乏 普 适 性 ， 导 致 对 广义 短路 比 


的 理解 相对 复杂 。 为 此 ， 本 文 利用 更 具 物 理 意义 的 


特征 子 系统 概念 


， 统 一 广义 短路 比 的 导出 原理 及 计 


算 方法 ， 主 要 讨 poner 
问题 1: 近似 原 大 规模 系统 模 态 的 低 维 解 直子 


系统 ( 即 特征 子 系统 ) 的 理论 基础 。 


问题 2， 基 于 特征 子 系统 ， 广 义 短路 比 及 其 临 


界 值 的 统一 导出 


原理 及 计算 方法 。 


3 ”特征 子 系统 及 其 近似 方法 


结合 定义 2 可 知 , 对 给 定 动态 系统 及 其 系统 矩 
阵 , 存在 低 维 特征 子 系统 与 广义 特征 和 矩阵。 其 中 ， 


特征 子 系统 可 近似 原 动 态 系 统 在 关注 模 态 下 的 稳 
定性 ， 且 维 数 更 低 ， 分 析 更 简单 。 


上 述 分 析 ) 


电路 的 关键 是 获得 特征 子 系统 , 即 求 


解 广义 特征 和 矩阵。 为 此 , 本 节 首 先 利用 矩阵 的 舒 尔 
补 变换 简化 原 系统 矩阵 ; 其 次 , 基于 引 理 1 获得 系 
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统 矩 阵 的 近似 特征 子 空间 , 并 得 到 对 应 的 广义 特征 
矩阵 及 低 维 特征 子 系统 。 
3.1 系统 矩阵 及 其 简化 

考虑 图 1 所 示 电力 系统 , 下 文 从 小 扰动 同步 视 
角 给 出 其 特征 子 系统 的 近似 求解 方法 ， 其 余 模 态 下 
的 特征 子 系统 求解 方法 类 似 。 

将 系统 看 成 由 两 部 分 组 成 , 一 部 分 为 等 效 网 络 ， 
包含 交流 网 络 与 m 台 构 网 装备 的 等 效 导 纳 ( 在 分 析 
频段 内 的 外 特性 由 戴 维 南 等 效 电路 主导 ， 其 等 效 时 
W Beg G=l,...m)); 男 一 部 分 包含 n 台 跟 网 装 


C 


备 动态 模型 。 式 (10) 所 示 的 ) 原 系统 矩阵 可 变换 为 ; 
TLGH) a 
HH: YO (s)=(B+Byp)@S(s) 
Vina (8) 
¥O()= Zou (8) 
Os mtr) 


Pue=| aem) || RO) FU) 
i 7(s) Ts) 
， Ber 为 m AIRE Ase BMT PH ZS R BCE FH BE 
的 修正 Tis). To(s). TOM Ta(s) SAA iE HE 
的 子 和 矩阵 ; Ib A 2 维 单位 矩阵 0， 0; 分 别 为 n、 
r 维 零 矩阵 。 
假设 未 接 装 备 时 电网 是 稳定 的 ， 则 根据 Schur 
ELI, det(Tne(s))=0 的 零点 由 det(Ti(s)- T2(s) Ty 
\(s)T3(s)) =0 决定 ， 故 det(T?,.(s))=0 决定 了 图 1 
所 示 系 统 的 稳定 性 ， 其 中 (推导 见 附 录 B): 


det (7, (s)-T, (s)T,;"(s)T;(s )) =der( T (s )) 
= det (YQ ()+¥2 (s))=0 


~ (s ) 


I 


(15) 


式 中 YR(s)= | 
Yer (8) 
Yu? (s)=Ba DG (s) 
B+B pr= K B, ‘| B', =(B,-B,B,'B, ) 
BB Pe 
H 中 ， B, ER", B ER”), BERO”, Ba 
ER DX) 。 
进一步 考虑 运行 点 的 影响 ， 将 剩余 台 跟 网 型 
装备 实际 工 况 与 其 动态 解 簿 2。 有 具体 地 ， 用 包含 装 
备 实际 工 况 信息 的 矩阵 (LCZ/PmQD AF Tres), 
TOjals) 可 转化 为 如 下 所 示 的 TOncds) 0 € M= 
TOrels) M 的 特征 函数 决定 了 图 1 系统 的 稳定 性 。 


第 刘 晨 三 等 ， 基 于 特征 子 
det(((C2/P (5)=aet(79 (s) 
oe et )(s))=0 
RP: ¥O(s)=(U es 
YR) =((U*/P) er, ia (S) = diae Yn) 
H, U?/P=diag(U?7/P;), Un Pi AAAS A i 的 端 


(16) 


ot 


0 


有 功 功 率 为 负 。 

JE 1: 当 装 备 在 额定 运行 点 , 其 端口 电压 为 1.0 
p.u.; 输 出 有 功 功率 等 于 其 额定 容量 Sn (二 1,...,n)， 
此 时 U?/P= diag(U2/P)-diag(Sw -Sw 故 额 定 运 
行 点 的 分 析 方 法 是 实际 工 况 下 的 特例 。 

3.2 ”特征 子 系统 近似 方法 

4 n 台 等 效 装备 动态 模型 同 构 时 ， 等 效 导 纳 矩 
阵 (U2/P)B’cn Qh KF AF BH 值 的 特征 子 空 
HAM 的 一 组 特征 子 空间 。 其 中 ， 当 系统 中 有 
(Gg<m 台 装备 有 功 功率 反 向 (为 负 信 ) 时 , (PPB sg 
OLA q 个 负 特 征 值 09， 将 特征 值 排序 为 1.4<...<4- 
OMH.. Ang M 为 最 小 正 特征 值 ，R(X,)、 
R(V,) G= -q…-1,1,.…1-9) 为 (U3/P)B rea Qo THE 
值 刀 对 应 的 特征 子 空 间 。 进 一 步 ， 根 据 定义 1 和 2 
并 类 比 式 (3)， 利 用 R(X) M R(Y,) (kE 人 9 
1,1,...,n-q}) 可 求解 原 系 统 的 多 个 特征 子 系统 ， 分 别 
近似 原 系 统 不 同 模 态 下 的 稳定 性 ， 总 结 为 如 下 引 理 
2， 其 具体 证 明 见 附录 C。 

引 理 2: 广义 特征 矩阵 二 满足: 

L= =Y" MX, =Y pry (5 )+4:¢(s) (17) 
RH: M=, QY pey (5)+(U°/P)B 86 (5) 
其 中 , Yrer DINAR E E ERRES; X= 
“Oh, YE ƏR, X Ve WAA a A FW 
(W7/P)B red@Qb RF AWA. ARES, xi yE 
R”, xe ye 分 别 为 (UWP)B rea KF ANA. EFF 
征 问 量 。 

当 装 备 动态 模型 存在 差异 时 , 构造 等 效 同 构 系 
统 , 其 系统 矩阵 的 特征 子 空间 可 近似 原 系统 矩阵 的 
特征 子 空间 ， 再 由 引 理 1 可 得 广义 特征 矩阵 五 的 
近似 。 具 体 地 ， 首 先 利 用 办 (XX ) AIRY, ) 构造 等 效 
同 构 系统 矩阵 Mi 为 : 

M, = 了 ww als )= Yor als )+F a. hom (s) (18) 
RP: Ya ton (8)=(U?/P) Bi, @E(s) 
了 ap 4 (s)=I, 8G, (s ) 


电压 、 输 出 有 功 功率 , 当 考 虑 装备 的 功率 反 向 时 ， 
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G, (s) = 了 | diag (WR Yn (s)) | XL 
=) Ps (U? [PY pr: (s), Pri = Xi Yri 
其 中 ， Thom (9) 为 对 应 本 的 等 效 同 构 系统 矩阵，Gx 
(5) 为 等 效 加 权 装 备 模 型 ;px 为 第 i 台 装 备 参 与 因子 ， 
Xkix Yri DBA Xk~x Vk 的 第 i~| 个 元 素 。 

其 次 ， 观 察 构 造 的 同 构 系 统 Mi 中 ， 装 备 动态 
是 原 系 统 中 所 有 接 入 跟 网 型 装备 动态 的 加 权 。 再 根 
据 引 理 1 和 2 TA, R(X) MRC) 仍然 是 Mi 的 
特征 子 空 间 ， 利 用 该 组 特征 子 空间 可 求解 履 的 广 
义 特征 矩阵 上 上, ， 表 示 为 : 

L, =Y" MX, =G,(s)+4,-¢(s) (19) 
A, -$(s)=¥i" ((U?/P)B,., 86 (s))X, 
局 (U°/P) Bax, -$(s) 

最 后 , 将 五 看 作 一 低 维 系统 的 系统 矩阵 , 那么 
该 低 维系 统 是 原 系统 的 一 个 特征 子 系统 , 可 近似 原 
系统 在 对 应 模 态 下 的 稳定 性 (由 于 保留 了 原 系 
统 和 矩阵 的 零 极 点 ， 故 还 能 反映 原 系统 的 鲁 棒 性 )。 
该 结论 总 结 为 定理 1， 具 体 证 明 见 附录 。 

定理 1: HR ASEM wm GER AGERE 
Mi 的 摄 动 ，Mi RETE ES or REEF Ze BIE 
似 相 等 ， 且 万 的 广义 特征 矩阵 五 可 由 式 (19) 计 算 。 
上 述 得 到 的 特征 子 系统 维 数 较 低 , 并 能 够 近似 
反映 原 系统 在 对 应 模 态 下 的 特性 。 因 此 ， 基 于 特征 
子 系统 简化 了 复杂 系统 的 分 析 。 

YE 2: 当 跟 网 装备 动态 模型 同 构 ， 等 效 同 构 系 
统 和 矩阵 与 原 系 统 和 矩阵 的 特征 子 空 s 间 完全 相同 ， 特征 
子 系统 稳定 性 与 原 系统 相同 ; 反之 ， 等 效 同 构 系统 
拖 阵 与 原 系统 矩阵 间 存 在 差异 , 广义 特征 矩阵 的 近 
似 与 真实 广义 特征 矩阵 存在 偏差 此 时 特征 子 系统 
近似 反映 了 系统 稳定 性 和 重 棒 性 。 

3.3 ”特征 子 系统 的 物理 意义 

据 上 节 可 知 , 利用 特征 子 空间 摄 动 可 求解 原 
系统 矩阵 的 广义 特征 矩阵 ， 对 应 的 特征 子 系统 可 
用 于 分 析 原 系统 稳定 性 ; 物理 上 看 , 特征 子 系统 是 


AF: 


一 个 包含 装备 与 网 络 动态 的 低 维 单 馈 入 系统 ， 本 
质 是 原 系统 在 模 态 坐标 系 下 的 一 个 解 耦 子 系统 。 特 
征 子 系统 的 物理 意义 及 求解 过 程 可 用 图 2 和 图 3 
总 结 。 


由 式 (19) 与 图 2 可 知 , 特征 子 系统 为 单机 无 穷 
大 系统 ， 由 跟 网 型 等 效 加 权 装 备 与 等 效 线路 构成 。 
其 中 ， 接 入 装备 的 特性 为 原 系统 中 跟 网 型 装备 的 加 
权 平 均 ， 反 映 系统 中 装备 动态 的 整体 特性 ;等 效 线 
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路 动态 计 及 了 构 网 装备 等 效 导 纳 等 因素 对 网 络 动 
态 的 影响 。 进 一 步 ， 从 传统 短路 比 看 ， 特 征 子 系统 
的 短路 比 取 值 为 对 应 等 效 线路 导 纳 ， 其 大 小 为 特征 
(E Ar 的 绝对 值 ， 可 用 于 表征 特征 子 系统 的 电网 强 
度 。 


实际 复杂 系统 


跟 网 型 装备 1 


5 
o C HNE 
pun YD 


” n+tm+r+l 外 部 电网 


特征 子 空间 | | 《下 稳定 性 
近似 。 ~ \Z 近 似 相同 


第 工 个 关注 模 态 下 的 特征 子 系统 


第 1 个 跟 网 型 
等 效 加 权 装 备 。 TRAM 


oti 


G,(s) 


SCR= |i; |=gSCR ~ 


第 个 关注 模 态 下 的 特征 子 系统 


第 n 个 跟 网 型 
等 效 加 权 装备 。 等 效 线路 


-E 


G.,(s) 


SCR = |2., | 


图 2 特征 子 系统 物理 意义 示意 图 
Fig.2 The schematic diagram of physical significance of 
eigen-subsystem. 
| 基于 特征 子 空间 摄 动 的 特征 子 系统 求解 
I 


构建 原 系统 矩阵 
Thes) 
见 式 (8)-(10) 


获得 原 系统 物理 模型 


和 矩阵 拆 分 、Schur 补 
得 到 简化 系统 矩阵 用 
见 式 (14)-(16) 


根据 稳定 模 态 ， 简 化 跟 
网 或 构 网 装备 动态 


基于 特征 子 
求解 特征 子 


空间 摄 动 
空间 近似 


近似 模 态 解 耦 


特征 子 系统 ， 解 耘 子 系统 
特征 子 系统 矩阵 L ( 见 式 (19)) 
注 : L, 为 广义 特征 矩阵 近似 


图 3 特征 子 系统 求解 方法 示意 
Fig.3 The schematic diagram of the calculation method of 
eigen-subsystem. 


还 值得 注意 的 是 ， 由 于 不 同 模 态 的 稳定 机 理 不 
同 ， 其 对 应 的 特征 子 空间 及 特征 子 系统 应 不 同 ， 对 
应 下 取 值 不 同 。 因 此 ， 实 际 问题 中 需 根据 具体 模 态 
及 稳定 机 理 ， 选 择 相应 的 特征 子 系统 ， 实 现 复杂 系 
统 的 简化 分 析 。 此 外 ， 由 于 本 文 聚焦 跟 网 型 装备 的 
主导 模 态 并 将 构 网 装备 等 效 为 对 地 支 路 ， 对 应 特征 
子 系统 仅 为 跟 网 型 装备 的 单机 无 穷 大 系统 ， 然 而 ， 
在 其 余 模 态 下 ， 对 应 特征 子 系统 形式 更 为 复杂 ， 包 
括 但 不 限于 单机 无 穷 大 系统 。 


4 基于 特征 子 系统 的 广义 短路 比 导出 原理 
及 强度 量化 方法 

本 节 提出 了 基于 特征 子 系统 的 广义 短路 比 统 
一 导出 原理 ， 并 给 出 了 相应 流程 和 几 种 典型 场景 下 
的 计算 方法 。 
4.1 特征 子 系统 中 的 电网 强度 和 广义 短路 比 指标 


根据 3.3 节 可 知 ， 特 征 子 系统 可 近似 表征 原 系 
统 的 稳定 性 , 并 由 等 效 加 权 装 备 动态 与 等 效 线路 动 
态 构成 ， 电 网 短路 比 大 小 为 特征 值 克 的 绝对 值 。 因 
此 ， 原 系统 强度 取决 于 特征 子 系统 动态 ， 当 等 效 加 
权 装 备 动态 一 定时 ， 原 系统 强度 取决 于 特征 子 系统 
的 短路 比 ， 即 等 效 线路 导 纳 。 

当 电网 电压 支撑 强度 低 时 ， 静 态 电压 稳定 及 由 
锁 相 环 带 来 的 小 扰动 失 稳 风险 较 高 ， 故 这 两 类 稳定 
问题 所 对 应 的 特征 子 系统 ( 即 所 1) 备 受 关注 。 具 
体 地 ， 在 该 模 态 下 ， 短 路 比 越 小 ， 系 统 失 稳 风险 越 
高 中; 那么 对 应 短路 比 为 如 (对 应 厂 1) 的 特征 子 系 统 
最 易 发 生 失 稳 , 该 短路 比 被 定义 为 “广义 短路 比 ”， 
如 式 (20) 所 示 。 因 此 ， 广 义 短 路 比 仅 是 基于 特征 二 
系统 的 一 个 具体 应 用 ， 其 余 重 要 的 特征 子 系统 ， 
能 够 反映 原 系统 在 其 余 模 态 下 的 特性 ， 而 且 都 能 
出 对 应 的 指标 。 然 而 ， 般 来 说 ， 基 于 短路 比 的 
析 方 法 适用 于 静 稳 及 小 扰动 同步 等 模 态 对 应 的 特 
征 子 系 统 ， 其 特点 是 这 类 问题 与 电网 电气 距离 强 相 
关 。 至 于 稳定 机 理 为 LC 谐振 的 高 频 振荡 对 应 的 特 
征 子 系统 ， 基 于 短路 比 的 分 析 方 法 可 能 不 再 适用 ， 
具体 见 文献 [7] 的 讨论 。 


gSCR=y;" (U?/P)By.,X,=Arin (C /P) B'a) 
-mo em 
=g" G JPY? Zu (U /有 


Ep, Ain (+) 表示 求 矩 阵 的 最 小 正 特 征 值 ，Zy 广 


= 


B red! o 
此 外 , 由 式 (20) 可 知 , 广义 短路 比 既 可 以 通过 和 拢 


当 闭 环 系统 的 特征 子 空间 与 开 环 系统 〈 即 仅 考虑 电 
Px) SRE RE) 的 特征 子 空间 相同 时 , 如 同 构 系统 中 ， 
两 者 结论 一 致 ， 反之 ， 基 于 特征 值 的 定义 仍然 是 一 
种 误差 较 小 且 便 于 实际 操作 的 计算 方法 。 
再 者 , 结合 式 (19) 和 (20) 还 可 知 ， 当 等 效 加 权 装 
备 / 场 站 动态 一 定时 , 根据 广义 短路 比 与 系统 强度 的 
正 相 关 性 , 能 够 确定 等 效 加 权 装 备 / 场 站 在 达到 临界 
稳定 状态 或 给 定性 能 要 求 所 需 的 最 小 短路 比 (简称 
“等 效 加 权 装 备 / 场 站 临界 短路 比 SCRo”)。SCRo 可 
反映 装备 / 场 站 动态 对 电网 强度 的 耐 受 能 力 , 恰恰 就 
是 广义 短路 比 临 界 值 CgSCRo。， 表 示 为 : 

SCR, = arg {det(G, (s,)+g5CR-¢(s,))=0} (21) 


其 中 ，arg{ } 表 示 求 方程 的 根 ， 54 为 等 效 加 权 装 备 在 
给 定性 能 要 求 下 解 得 系统 主导 特征 根 。 

综 上 ， 复 杂 系 统 的 广义 短路 比 统一 导出 过 程 为 : 
1) 利用 等 效 同 构 系 统 特征 子 空间 ， 求 解 原 系统 矩 
阵 近 似 的 广义 特征 矩阵 , 获得 能 够 近似 原 系 统 动态 
特性 的 特征 子 系统 ; 2) 基于 特定 的 特征 子 系统 , 计 
算 其 短路 比 指标 ( 即 广义 短路 比 ) 及 装备 / 场 站 临界 
值 ( 等 于 广义 短路 比 临界 值 ), 进而 形成 融合 电网 ) 
义 短 路 比 与 装备 / 场 站 临界 短路 比 的 量化 方法 。 其 中 ， 
广义 短路 比 仅 与 等 效 网 络 动态 和 接 入 装备 的 容量 
有 关 、 临 界 值 仅 与 等 效 加 权 装 备 / 场 站 的 动态 有 关 ， 
故 广义 短路 比方 法 保留 了 传统 短路 比方 法 所 具有 
的 网 络 动态 与 装备 动态 分 离 量化 的 优势 。 
42 ”典型 场景 下 广义 短路 比 计算 方法 

结合 3.1 节 中 系统 矩阵 的 简化 方法 、 公 式 ( 见 式 
(15)-(16)) 以 及 4.1 节 中 广义 短路 比 的 计算 公式 ( 见 
(20)), 本 节 给 出 三 种 典型 场景 中 系统 矩阵 的 具体 形 
式 ， 由 此 实现 广义 短路 比 的 计算 。 
4.2.1 考虑 装备 在 非 额定 运行 点 的 广义 运行 短路 比 

考虑 图 1 所 示 系 统 仅 包 含 n 台 采 用 不 同 控制 策 
略 、 参 数 的 新 能 源 变 流 器 ，m=0。 其 中 ， 新 能 源 变 
流 器 在 不 同 的 非 额定 运行 点 ， 输 出 有 功 功率 及 端口 
电压 分 别 为 0.2~1p.u. 和 0.9~1.1p.u. 不 等 231]。 根 据 式 
(15)-(16), (20), 即 可 求 得 该 场景 下 简化 后 的 系统 矩 
阵 及 广义 短路 比 。 具 体 过 程 为 : 

首先 ， 考 虑 无 构 网 装备 等 效 导 纳 对 网 络 动 态 的 
修正 , 即 Bzsr=0n+, 将 Bisr=0n4r 代 入 式 (15) 并 令 m=0， 
可 求 得 仅 包 含 h 台 装备 接 入 母线 的 网 络 等 效 导 纳 矩 
阵 Breas 其 次 ,考虑 将 nn 台 装 备 的 实际 工 况 与 其 动 
SRO YO (s) 5 U?/P 如 式 (22) 所 示 ; w, K 


阵 特 征 值 定义 ， 也 能 通过 闭环 系统 的 特征 向 量 得 到 。 
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Breas YQ (s) 与 02/ 忆 分 别 代 入 式 (16)、(20)， 可 求 
得 该 场景 下 的 系统 矩阵 与 广义 短路 比 。 

YM (s) = diag (¥isc;(8)),U?/P = diag (U? /P) (22) 
其 中 ，Yysci(s) 为 第 i 台新 能 源 变 流 器 导 纳 传递 函数 
FERE, 二 1,.….,n; Un Pi 分 别 为 装备 i 的 端口 电压 、 
输出 有 功 功 率 。 
4.2.2 考虑 装备 有 功 功率 反 向 的 广义 短路 比 

根据 装备 的 物理 特性 与 系统 需求 ， 部 分 装备 可 
反 向 从 电网 吸收 功率 ， 如 储 能 系统 与 柔性 直流 输电 
等 。 不 失 一 般 性 地 ， 以 常规 储 能 系统 为 例 ， 考 虑 图 
1 所 示 系 统 包含 n 台 采 用 不 同 控制 策略 、 参 数 的 新 
能 源 变 流 器 ，m 台 常 规 储 能 系统 。 其 中 ， 装 备 均 运 
行 在 额定 工 况 ; 新 能 源 变 流 器 主要 向 电网 发 出 有 功 
功率 ， 储 能 变 流 器 主要 向 电网 发 出 (放电 ) 或 吸收 有 
Th (充电 )。 根据 式 (15)-(16)、(20), 即 可 求 得 该 场景 
下 简化 后 的 系统 矩阵 及 广义 短路 比 。 具 体 过 程 为 : 
首先 ， 考 虑 无 构 网 装备 等 效 导 纳 对 网 络 动 态 的 
修正 ， 即 Bisr=0w4m+:， 将 Bipr=On+mir 代入 式 (15)， 
可 求 得 仅 包 含 ntm 台 装 备 接 入 母线 的 网 络 等 效 导 
内 矩阵 Breas HK, SUR ntm 台 装 备 的 额定 容 
量 与 其 动态 解 而 ， 其 中 ， 当 储 能 变 流 器 充电 ， 其 额 
定 容量 需 加 上 负 号 , 例如 ， 考虑 m 台 储 能 均 工作 在 
RERS, YO (s) 与 UV/P 如 式 (23) 所 示 ; 最 终 , 将 
B’reas Yigg (s) 与 D2/ 忆 分 别 代 入 式 (16)、(20)， 可 求 
得 该 场景 下 的 系统 矩阵 与 广义 短路 比 。 

ye (s) 上 ba (Ysa (s)) 


五 


ALS 


tin 


diag (Veg (s j 


diag (Sxt) 
-diag( Sxia) 


其 中 ，Ypcs(s) 为 第 j 台 储 能 变 流 器 的 导 纳 传递 函数 
HERE, j=1,... m; Snin Snip PRAHA i 台新 能 源 
变 流 器 及 第 j 台 储 能 变 流 器 的 额定 容量 。 
4.2.3 考虑 构 网 型 装备 影响 的 广义 短路 比 

考虑 图 1 所 示 系 统 包含 n 台 采 用 不 同 控制 策 
略 、 参 数 的 新 能 源 变 流 器 ，zz 台 构 网 型 储 能 系统 或 
调 相 机 。 其 中 ， 装 备 均 运行 在 额定 工 况 ; 调 相 机 和 
外 特性 近似 为 电压 源 且 惯 量 系数 较 大 的 构 网 型 变 
流 器 能 够 为 系统 提供 电压 支撑 能 力 ， 在 锁 相 环 主导 
的 小 扰动 同步 稳定 问题 中 能 够 等 效 为 含 内 电抗 的 
电压 源 。 根据 式 (15)-(16)、(20)， 即 可 求 得 该 场景 下 
简化 后 的 系统 矩阵 及 广义 短路 比 。 有 具体 过 程 为 : 
首先 , 考虑 m 台 构 网 装备 等 效 导 纳 对 网 络 动态 
的 修正 ， 即 Brer= diag(0 diag(Ber) 0 ; )， 将 其 代入 
式 (13)， 可 求 得 仅 包含 n 台 装 备 接 入 母线 的 网 络 等 
效 导 纳 和 矩阵 B”on; 其 次 ,考虑 将 n 台 装 备 的 额定 容 


(23) 


2/Pp_S- 一 
U [P = Sw = 
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ESRAS, YO (s) 与 V0?/P 如 式 (24) 所 示 ; 最 
终 ,将 Breas YO (s) 5 U2/P A BARA SK(16), (20), 


可 求 得 该 场景 下 的 系统 矩阵 与 广义 短路 比 。 
YẸ (s )= diag( Yale )) 
U’/P =S} =diag (S; ') 


Co, Bey 为 第 j 台 构 网 型 变 流 器 或 调 相 机 的 等 效 
导 纳 ，J=1,...,m，Blsr=diag(SvswBp 或 1/xa;”)> Xaj” 


(24) 


Ni 


工程 学 R 


为 第 j 台 调 相机 的 等 效 内 电抗 ， 即 次 暂 态 电抗 ， 
1/SvswBp 为 第 j 台 构 网 型 变 流 器 的 等 效 内 电抗 ， 
Syy 为 第 j 台 构 网 型 变 流 器 的 额定 容量 ，Bp 为 单位 


容量 下 构 网 型 变 流 器 的 等 效 导 纳 。 


为 直观 对 比 上 述 3 种 典型 场景 下 原 系统 矩阵 的 


构成 ， 表 1 总 结 了 上 述 公 式 。 此 外 ， 


上 述 场 景 对 系统 强度 的 作用 机 理 ， 
参考 文献 [9][10][11][19][20]。 


表 1 典型 场景 下 原 系统 矩阵 构成 及 其 对 系统 强度 的 作用 机 理 


Tab.1 The elements required for calculation of system strength and their physical significance under typical 


#1 还 给 出 了 


详细 讨论 可 分 别 


scenarios 
ARSE M=YO (s)+¥O (s) 
7 )(s)= ((U?/P)®1, ) Y (s) YË (s)=(U?/P) Ba @¢(s) 
典型 场景 影响 机 
= Bran BiB ENE YOO AA ai 
Yeas) UP 等 效 导 纳 矩阵 
B 构 网 装备 等 效 导 纳 Bysr 
非 额 定 运行 点 diag(Yysci(s)) diag(U7/P) | 包含 所 有 装备 Ors» 
有 功 功 ZA 癌 diag( Yrsci(s) Y, pcs(S)) Syintmy! 及 中 间 节 H? 点 的 Ormrr 改变 gSCR 或 SCRo 
构 网 型 装备 diag(Yysci(s)) Syn" 网 络 导 纳 矩阵 diag(0, diag(Bran) 0, ) 
压 稳 定 外 的 其 它 稳定 问题 。 


43 ”系统 强度 量化 步骤 

为 方便 实际 应 用 ， 本 节 给 出 基于 广义 短路 比 的 
系统 强度 量化 步 又， 具体 如 下 
) 根据 装备 物理 特性 将 装备 动态 分 别 看 作 跟 
网 型 装备 动态 或 对 网 络 动 态 的 修正 ; 

2) 根据 网 络 参数 与 装备 特点 ， 基 于 式 (15) 得 到 
age rap alee ia et 
Bier FEF AUDIRE OCS FE ME (O?/P) Breas 

i abe ERE 
特征 值 或 奇异 值 ， 即 广义 短路 比 ; 

4) 根据 式 (21) 或 (25) 求 解 SCRo。 有 具体 地 ， 前 者 
可 通过 解析 计算 或 仿真 测试 等 方式 获得 门 ， 得 到 的 
SCRo 精度 较 高 ; 后 者 则 需 结合 实际 工程 经 验 , 利用 
参与 因子 ( 见 式 (18)) 计 算 获 得 。 


SCR, = piSCR, ， (25) 


式 中 ，SCRo ;为 第 i 人 台 装 备 达 到 临界 稳定 所 需 的 最 
小 短路 比 。 

5) 根据 式 (12) 求 解 6%， 量 化 强度 。 

值得 注意 的 是 ， 当 系统 为 黑箱 模型 时 ， 等 效 网 
络 动 态 与 实际 工 况 矩阵 可 通过 扰动 辨识 方法 人 2 与 
PMU 量 测 等 方法 获得 数据 。 

综 上 所 述 , 通过 引入 特征 子 系统 , 统一 了 多 种 
场景 下 广义 短路 比 及 临界 值 的 导出 原理 及 计算 方 
法 ， 具 有 优势 为 : 1) 保留 物理 特性 。 特 征 子 系统 
保留 了 装备 与 网 络 的 动态 特性 ; 2 ) 方 法 通用 且 拓 展 
性 好 。 不 仅 适 用 于 异 构 多 馈 入 系统 的 复杂 场景 ， 而 
且 有 潜力 推广 应 用 到 除 小 扰动 同步 稳定 和 静态 电 


一 


) 


5 算 例 分 析 及 验证 


首先 基于 4 人 馈 入 系统 验证 特征 子 系统 求解 方法 


TT 


及 其 在 描述 原 系统 稳定 性 的 有 效 性 ; 


HK, AAAS 


实际 新 能 源 基地 的 电磁 暂 态 模型 ， 验 证 广义 短路 比 
指标 在 新 型 电力 系统 中 的 有 效 性 ， 并 分 析 装 备 在 非 
额定 运行 点 、 有 功 功率 反 向 和 包含 构 网 型 装备 对 指 


标 计算 与 系统 强度 的 影响 规律 。 


5.1 ”特征 子 系统 及 其 求解 方法 的 有 效 性 

为 验证 特征 子 系统 在 描述 原 系 统 稳 定性 的 有 
效 性 ， 搭 建 4 馈 入 系统 电磁 暂 态 仿真 模型 拓扑 如 
图 4 所 示 。 其 中 节点 1、4 为 常规 、 构 网 型 储 能 系 
统 接 入 节点 , 节点 2、3 为 新 能 源 变 流 器 接 入 节点 ， 


J 


图 44 馈 入 系统 


装备 控制 策略 、 参 数 及 网 络 参数 如 附 表 El1-E3 所 示 。 


Fig.4 The topology of a four-converter system 


LAS, FS RAE EL PLS A I 
路 电感 ， 比 例 参数 k=l 1.5. 2, 3, 


E3 所 示 所 有 线 
4， 其 中 ， 常 规 


储 能 变 流 器 从 电网 吸收 功率 ， 构 网 3 


型 变 流 器 向 电网 


发 出 功率 , 装备 实际 工 况 如 表 2 所 示 。 以 此 为 基础 ， 
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首先 ， 分 别 基 于 特征 子 空 间 摄 动 理论 求解 特征 子 系 
统 ; 其 次 ,对比 原始 4 馈 入 系统 和 特征 子 系统 的 主 
导 特 征 根 ， 验 证 特征 子 系统 的 有 效 性 。 

首先 ， 以 算 例 1 为 例 ， 给 出 特征 子 系统 的 具体 
求解 过 程 。 有 具体 步骤 为 : 构建 系统 矩阵 ， 考 虚构 网 
型 变 流 器 对 网 络 的 修正 ， 将 系统 矩阵 降 阶 ， 并 用 


(UIP) Qh EREE: HARARE: 
det (T,) (s)) =det(¥in(s)+¥L)(s))=0 26) 
-55 14 14] 
式 中 : (U?/P)B,=|2 7 2 
4 4 16| 


' fr U2 U? 
YO (9) =a P Yocs (s) p iso (s) pe 加 
APF, Frsci (8)、Yyscz (3)、Ypcs(s) 分 别 为 新 能 源 变 
流 器 与 常规 储 能 变 流 器 的 导 纳 传递 函数 矩阵 ，ZP/P 
=diag(U1?/P U2?/P2 L32WP3)。 
12=17.10, A1=-53.62, EF Z R(X). R(Y,) g= 
-1,1,2) 满 足下 式 , HF X=x8h, Yay h, x, 
yi DNA IP) Brea KF yA ERE E. 
Y} ((U°/P)B a @1,)X=4,@L, 27) 
it ASE RERA ZAST REE E, BEZA 
特征 子 系统 的 系统 矩阵 ， 表 示 为 : 
L SYST (s)X;=G,(s)+4;-¢(s) (28) 
式 中 : Pi = Xe 


— 


2 2 


U. U2 U. 
G, = Pi pS (s)+ Dj, pra (s)+p; p Ysc (s) 


3 
HP, 1,2,3, xj. vi PRA xj yE i 个 元 素 。 
表 2 算 例 的 比例 参数 及 实际 工 况 
Tab.2 The proportional parameters and actual op- 
erating conditions of examples 


> 


算 例 | k U P gSCR 
1 1 | 1.05, 0.99, 1.03 5.19 
2 1.5 1.04, 0.97, 1 3.39 
3 2 | 1.03, 0.95, 0.99 | -0.2、 1, 0.6 2.47 
4 3 | 1.01, 0.95, 0.98 1.67 
5 4 | 0.99, 0.93, 0.97 1.21 


图 5 系统 等 效 前 后 主导 特征 根 对 比 图 


Fig.S Comparison of dominant roots before and after 


equivalence 

进一步 ， 结 合 4.1 节 分 析 可 知 ， 理 论 上 最 小 正 
特征 值 Ar 对 应 的 特征 子 系统 能 够 近似 表征 原 系统 
的 小 扰动 同步 特性 。 为 验证 该 结论 , 图 5 给 出 4 馈 
入 系统 和 1 对 应 特征 子 系统 的 主导 特征 根 。 可 以 看 
出 ， 随 着 电网 强度 的 变化 ， 原 始 系统 和 该 特征 子 系 
统 的 特征 根基 本 一 致 ， 验 证 了 该 特征 子 系统 可 近似 
原 系 统 的 小 扰动 同步 特性 。 
5.2 ”基于 广义 短路 比 的 系统 强度 量化 

搭建 菜 实际 新 能 源 基 地 的 电磁 暂 态 仿真 模型 ， 
如 图 6 所 示 。 该 系统 共 包 括 4 个 区 域 共 27 个 风力 
发 电机 组 、2 个 常规 储 能 系统 和 2 个 构 网 型 储 能 系 
统 ， 其 中 风力 发 电机 组 或 储 能 系统 的 下 标 数字 代表 
其 内 部 采用 变 流 器 编号 ; 图 中 节点 1~54 表示 装备 
接 入 节点 , 其 中 节点 43, 54 为 构 网 型 变 流 器 接 入 节 
点 ， 节 点 14、15、28、29 为 常规 储 能 变 流 器 接 入 节 
， 其 余 节 点 为 新 能 源 变 流 器 接 入 节点 ; 节点 
55~92 表示 中 间 节 点 ; 节点 93 为 场 站 外 电网 等 值 
节点 。 有 具体 地 ， 考 虑 如 下 四 个 算 例 情况 ， 算 例 的 装 
备 控制 策略 、 参 数 及 网 络 参 数 均一 致 ， 如 附 表 El- 
E3 所 示 ; 装备 实际 工 况 如 附 表 E4 所 示 。 

算 例 1: 系统 中 装备 均 运行 在 实际 工 况 ， 即 输 
出 有 功 功率 及 端口 电压 往往 偏离 额定 值 ， 其 中 ， 常 
见 储 能 从 电网 吸收 功率 ， 其 余 装 备 回 电网 发 出 功率 ; 
算 例 2: 在 算 例 1 的 基础 上 将 装备 工 况 均 调整 
为 额定 工 况 ， 其 中 ， 常 规 储 能 仍 从 电网 吸收 功率 ， 
其 余 装 备 癌 电网 发 出 功率 ; 
算 例 3: 在 算 例 1 的 基础 上 将 常规 储 能 变 流 器 
的 工作 模式 改变 为 发 出 功率 ; 

算 例 4: 在 算 例 1 的 基础 上 切除 构 网 型 变 流 器 ， 
即 系统 中 构 网 型 变 流 器 接 入 节点 悬空 。 


© i © | | ® | 
p ao i 2) | on Ed 
= a Bh ats osm 


a 


v 


t | 
se 39 
m2. 
i r 220kV 
y 口 a 0.69 kV 
p Meos Liners.s0 = 风电 机 组 
Lio Eo 常规 储 能 
Lines1,78 
42 5 vx 
ote 构 网 型 储 能 
Lines2,78 


图 6 某 实 际 新 能 源 基地 拓扑 图 
Fig.6 The topology of a practical renewable energy plant 


首先 ， 按照 上 述 指标 计算 流程 分 别 解析 计算 四 
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个 算 例 中 系统 的 广义 短路 比 ， 其 次 ， 基 于 式 (21) 并 


中 


国 电 


机 


通过 搭建 等 效 加 权 单 机 仿真 或 解析 方法 得 到 四 个 


算 
gSCR 45 SCRo Ft A Ata xe 
计算 po%， 以 评估 系统 强度 。 


侈 中 系统 的 SCRo; 最 后 ， 基 于 式 (11)， 通 过 比较 
FE, 进一步 ,基于 式 (12) 


表 3 给 出 评估 结果 。 可 以 看 出 ， 算 例 1、3、4 


中 gSCR 均 大 于 SCRo， 其 中 ， 


A 


例 1 “A B%>20%, 


满足 小 扰动 稳定 裕 度 需求 ; SF 


例 2 中 gSCR=2.39 


<2.63=SCRo, 系统 失 稳 。 该 判 稳 结果 与 图 7(a)(b) 所 


示 算 例 1、2 中 装备 端口 电压 的 时 域 仿真 趋势 一 致 ， 
验证 了 基于 广义 短路 比 的 量化 方法 的 有 效 性 。 
表 3 算 例 的 指标 数值 
Tab.3 The values of indexes under four scenarios 
算 例 gSCR SCRo p% 
1 3.28 2.67 22.8% 
2 2.39 2.63 -9.1% 
3 3.03 2.68 13.1% 
4 3.02 2.69 12.3% 


gSCR=3.28>2.67=SCRo 


P i 
& 1.05 
= 
1 
0.95 — 
0 02 04 06 0.8 1 1.2 
t/s 
(a) 
12l gSCR=2.39<2.63=SCR 
A 
a 1 
D 
0.8 
0.6 moa 4 aa 
0 02 04 06 08 1 1.2 
t/s 
(b) 
L2f <a 
gSCR=2.67=SCRo 
= 1.1 
fay 
S į 
0.9 
0 0.2 04 06 08 1 1:2 
t/s 
(c) 
图 7 装备 端口 电压 的 时 域 波 形 


Fig.7 Time domain waveform of terminal voltages of 
devices in a practical renewable energy plant under differ- 
ent Operating conditions 


此 外 ,图 7(c) 给 出 算 例 1 中 系统 在 gSCR= SCRo 
时 装备 端口 电压 的 等 幅 振荡 波形 ， 验 证 了 gSCR= 
SCRo 在 描述 系统 稳定 边界 的 有 效 性 , 此 时 系统 网 络 
参数 如 附 表 E3 所 示 。 


ay 


C 程 qk 


braa 


5.3 ”广义 短路 比 的 影响 因素 
5.3.1 ， 非 额定 运行 点 的 影响 

考虑 对 比分 析 装 备 在 实际 工 况 ( 算 例 1) 与 额定 
工 况 的 场景 ( 算 例 2)。 
首先 , 在 指标 计算 过 程 中 , 算 例 1、2 的 区 别 仅 
FEF IESE UVP 的 不 同 。 具体 为 : 额定 工 况 下 UVP= 
S Ns21，SN5 为 52X52 AXA FORE, TOR AE 
定 容量 ， 当 装备 吸收 功率 ， 额 定 容 量 需 加 上 负 号 ; 
实际 工 况 下 LP/P= diag(U2/Pi)ss,，diag(U2/Pis; 为 52 
x 52 WY FA FER. 

其 次 , 如 表 3 所 示 , 算 例 1、2 中 gSCR 及 SCRo 
都 随 着 工 况 而 发 生 了 变化 ， 说 明了 实际 工 况 对 系统 
强度 的 影响 。 
5.3.2 ”有 功 功率 反 向 的 影响 

考虑 对 比分 析 常 规 储 能 变 流 器 吸收 有 功 功 率 
( 算 例 1) 与 发 出 有 功 功率 的 场景 ( 算 例 3)。 
首先 ， 根 据 常 规 储 能 变 流 器 在 充电 或 放电 工作 
模式 下 的 物理 特性 可 知 ， 在 指标 计算 过 程 中 ， 算 例 
1、3 的 主要 区 别 在 于 和 矩阵 UP 的 不 同 。 有 具体 为 ， 
常规 储 能 变 流 器 吸收 功率 时 , U2/P 中 与 其 接 入 节点 
对 应 的 对 角 元 为 负数 ， 使 得 等 效 导 纳 和 矩阵 有 对 应 的 
负 特 征 值 ， 而 发 出 功率 时 ， 相 应 对 角 元 为 正 数 ， 对 
应 特征 值 均 为 正 数 。 

其 次 , 比较 两 算 例 中 指标 及 临界 值 。 可 以 看 出 ， 
相 较 于 储 能 变 流 器 工作 于 放电 模式 ,功率 反 向 后 系 
统 的 gSCR 得 到 一 定 提升 ， SCRo 稍微 改变 , 说 明 储 
能 变 流 器 工作 于 充电 模式 会 提升 广义 短路 比 并 改 
AE SCR。。 
5.3.3 构 网 型 装备 的 影响 

考虑 构 网 型 变 流 器 接 入 ( 算 例 1) 与 不 接 入 的 场 
景 ( 算 例 4)。 
先 ， 根 据 构 网 型 装备 的 物理 特性 可 知 ， 在 指 
标 计算 过 程 中 , 算 例 1、4 的 主要 区 别 在 于 装备 对 等 
效 网 络 动态 的 修正 Bi 不 同 。 有 具体 为 ， 考 虑 构 网 型 
变 流 器 接 入 时 , Bysr 包含 构 网 型 变 流 器 的 等 效 导 纳 ， 
会 修正 网 络 动态 ， 不 接 入 时 Ber 中 元 素 均 为 0。 

其 次 ， 比 较 算 例 1、4 中 广义 短路 比 及 SCRo. 
可 以 看 出 ， 接 入 构 网 型 变 流 器 后 系统 广义 短路 比 得 
到 一 定 提升 ， SCRo 稍微 改变 , 说 明 构 网 型 变 流 器 的 
接 入 会 提升 广义 短路 比 并 改变 SCR。。 
6 结论 与 展望 

本 文 给 出 了 多 馈 入 系统 的 特征 子 系统 定义 、 物 
理 意义 及 计算 方法 ， 并 提出 复杂 场景 下 广义 短路 比 
及 其 临界 值 的 通用 导出 原理 与 计算 方法 ,包含 在 非 
额定 运行 点 、 有 功 功率 反 向 和 包含 构 网 型 装备 等 多 
种 场景 。 特 征 子 系统 保留 了 交流 网 络 与 跟 网 型 装备 


I 


fF 


的 动态 特性 ， 存 在 与 之 对 应 的 低 维 物理 系统 ， 故 由 
此 导出 的 广义 短路 比 计 算 和 系统 强度 量化 方法 具 
有 较 严 说 的 理论 基础 ， 也 有 较 清 晰 的 物理 意义 。 

本 文 聚焦 跟 网 型 装备 主导 的 稳定 模 态 ， 其 特征 
子 系统 为 跟 网 型 装备 的 单机 无 穷 大 系统 ; 其它 模 态 
下 ， 如 构 网 型 装备 控制 特性 不 能 忽略 时 ， 特 征 子 系 
统 还 需要 保留 构 网 型 装备 ， 其 特性 需 深入 研究 。 进 
一 步 ， 特 征 子 系统 可 能 是 电力 系统 的 基本 单元 , 包 
含 但 不 限于 单机 无 穷 大 系统 ， 其 基本 理论 和 方法 后 
续 还 需 进 
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附录 A 引 理 1 的 举例 
以 两 维和 矩阵 为 例 说 明 引 理 1。 设 M 与 详 分 别 


为 : 
1 2 1 2 | 
-| erm 7 errdiog(s) (29) 
其 中 ，gi;E R， 数 值 较 小 ，i=1,2,3,4。 
根据 式 (29), 利用 M 的 两 个 特征 子 空间 可 求解 
应 广义 特征 矩阵 的 近似 ; 
L ~Y" MX, 
0.756, +0.24¢. (30) 
=~0.37@1 ! 
al | 
L, =Y; MX, 
.246,+0.7 31 
Soba ANPR ea 
0.248, +0.768, 


ai ee 057] @1,,X,=[-042 -0.91] 81, 
4 =[.0.92 0.42]@1,, ¥/ =[-0.57 -0.84]®L, 


其 中 ， ee see 重 特征 值 -0.37 
55.37 的 右 特征 向 量 , Yi. Yo PAI ele At MAE 


特征 向 量 ， ican E 子 空间 。 
如 果 er. en 63, 64 DHA 0.4, 0.3. 0.2. 0.1, 


工 程 学 报 
必 广 义 特征 矩阵 及 其 近似 结果 分 别 为 
-0.021 5.621 
5| ya A oe 
-0.022 5.618 
=| | | ce 
附录 B 式 (15) 的 详细 推导 
原 系统 的 系统 矩阵 为 : 
oral T) TCs) 
ATO) nO icy 


stein, 1 (s)=¥2(s)+B, 8¢(s).1,(5) =B, 86s) 
1,(s)=B, 86(8)sT,(s)=B, ®6(s) 
FFE RE RIA RADE, m: 
T (s)-7,(s)T- (s)T(s) 


35 
=¥(2(s)+Bi, ec 人 pa 

附录 C 51 2 的 证 明 
根据 特征 值 与 特征 子 空间 定义 可 知 ， 
(ZLP2/ 门 Byez QLA n ZERIE y IAA 
特征 子 空间 R(X )、 R(Y,), 中 ， Xj =x; ® h, 
Y= yi" @ Ly Xn YW A(O/P)B rea 8h 
关于 AMA. AMER. Al, MES 1 式 


(3), FA R(X,) AUR(Y,) RE WP)B ea Ob ih n 
个 广义 特征 矩阵 ; 
L=Y"((V’/P)B', ®L)X,=% 1 (36) 
进一步 ， 当 原 系统 矩阵 让 中 等 效 装备 模型 相似 ， 
My XA Yi eo pe: 
MX, =X,L,,Y"'M=LY? 
式 中 : L= Yora (S )+4, g(s ) 
比较 式 (1)、(37) 并 结合 特征 子 空间 定义 (定义 1) 
ays, R(X). R(X) 和 工分 别 为 原 系 统 矩 阵 的 特 
征 子 空间 及 广义 特征 和 矩阵, 因此 , 令 X=(Xq... Xna) 
Y=(¥Vig... Yn) 利用 XX 了 可 将 原 系 统 矩 阵 转换 为 : 
L=Y" MX =diag(L,) (38) 
EX All Y AE BEIJAR PEJE, ie det ( M )=0 
5 de(L)=0 %1, WG det(L, )=0 等 价 。 因此 ， 对 


应 的 特征 子 系统 (五 为 特征 子 系统 矩阵 ) 可 近似 原 
系统 的 稳定 性 。 证 毕 。 
附录 D 定理 1 的 证 明 

证 明 : 利用 R(X,) RO) 可 求解 等 效 同 构 系 
GERE Mt 的 广义 特征 矩阵， 满足 : 


(37) 


基于 特征 子 系统 的 广义 短路 比 导 上 
L, =Y"M,X, =G, (s) +A, -¢(s) 


(39) 


co 
z 
H 


里 及 计算 方法 13 
RE 仿真 电网 参数 


Table E3 Parameters of Simulations 


Fee Py y be / E a7 
结合 式 (39)、(19), SSC ARACHIS IE RE URAK FM SMa) 
y ck. + fy = = = = = 
阵 是 原 系统 广义 特征 矩阵 的 近似 。 根 据 引 理 1， 利 Xıs5=0.20 AeA” X30.60 Xiz040  Xz3=0.40 
P AEZ 区 X13=0.40 X3,4=0.20 
H + y 3H = 
I R(X, ) ALR, ) 可 得 广义 特征 矩阵 的 近似 * 证 毕 。 Sie IE EA SROUO pay 
附录 E X55, 86=0.24 X56.86=1.00 X57.86=2.40 X58.86=3.30 X59,86=2.90 
X60,86=2.00 X61,86=1.70 X62,86=0.03 X63,87=1.10 X64.87=3.10 
KEL 新 能 源 及 储 能 变 流 器 控制 参数 eT ATE 
Table E1 Control parameters of grid-following converters X70,87-0.02 — X71,8=2.50 X72882.20  Xnss=1.50  X74,88=1.30 
内 容 新 能 源 变 流 器 储 能 变 流 器 X75,88=0.81 X76.88=1.30 X77,38=0.81 X78,90=0.19 X79,90=2.20 
Xgo,g9=0.19 Xs1,89=2.20 Xg2,91=1.50 Xg3,91=1.30 Xs4,91=1.50 
Ax ps ZA bos HT az oe 
ane AMEE 1500kVA, 620V Xs5s,91=1.30  Xs6,87=1.00  Xs7,92=3.10  Xg8,92=4.40  Xs9,90=7.00 
E HOS Xo0,92=7.80  Xo1,92=0.87  Xo2,93=6.90 变压器 电抗 : 64.0 
en er 0.68MVA(14,15) pe ne ek £ 
装备 额定 容量 1.5MVA 算 例 1 在 gSCR=SCRo 时 的 电网 所 有 参数 是 上 述 电网 参数 的 1.1 倍 
0.72MVA(28,29) 
直流 侧 电 容 、 电 压 0.1089pu,LIEV  0.1089p.u., 1.5kV 表 E4 某 实 际 新 能 源 基 地 4 个 算 例 中 的 装备 实际 工 况 
外 环 PI 0.5,40( 直 流 电 压 ) 0.5,20( 功 率 ) Table E4 The operating conditions of devices in four examples 
滤波 电感 、 电 容 0.05、0.0006 0.05、0.0006 of Practical Renewable Energy Plant 
ss ek 0.5,15(14,15) 算 例 运行 工 况 
we ; 
电流 内 环 PI 00 0.3,10(28,29) F 1-13:0.7;16-27:0.8;30-38:0.6;39-42:0.75; 44-53:0.9; 14-15: - 
0.45;28-29: -0.48;43,54:1.00; 
paa RB, > st, ; 
| 局 入 系统 : ogee I y | 1:13, 16-27, 30-41, 44-48 :1.01;14-15, 28-29:1.03;42,43,54:1.00; 
某 实际 新 能 源 场 站 : 49-53:1.02 
0,10500(1-13); 21,10300 P. 1.00 
锁 相 环 PI a 40,24000 2 leg 1.00 
(16-27); 26,10400 (30-38); —— 
‘ P 14-15: 0.45;28-29: 0.48; 其 余 见 算 例 1 
6: 102003354244548); U 14-15, 28-29:1.01， 其 余 见 算 例 1 
31,10200(49-53) P 见 算 例 1 
电压 前 馈 时 间 常 数 0.00001 t Fy 见 算 例 1 
表 E2 构 网 型 变 流 器 控制 参数 seer 
Table E3 Control parameters of grid-forming converter 收 稿 日 期 : 
内 容 数值 作者 简介 : 


装备 额定 容量 


滤波 电感 、 滤 波 电容 、 电 阻 
电压 外 环 、 电 流 内 环 PI 参数 


构 


= 
as 


电压 前 馈 滤波 时 间 常 数 


1.5MVA 


系数 、 阻 尼 系 数 2,20 


0.001 


0.05、0.0005、0.02 
0.8,12; 6,10 
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Derivation Principle and Calculation Method of Generalized Short-circuit Ratio Based on 


Eigen-Subsystems 


XIN Huanhai!, LIU Chenxi!, HUANG Linbin!, GAO Huisheng!, JU Ping! 
(1. College of Electrical Engineering, Zhejiang University) 
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With the large-scale integration of power electronic 
devices, voltage support strength(referred to as “system 
strength”) decreases, thus threatening the system's safe 
and stable operation. In both the homogeneous and heter- 
ogeneous scenario (i.e., power electronic devices inte- 
grated into the system have the same dynamics or weakly 
different dynamics), system strength can be quantified by 
the gSCR and SCRo. However, there is a lack of unified 
derivation principles and methods. To this end, this paper 
unified these based on the concept of eigen-subsystems. 

Specifically, we consider a power system integrated 
with types of devices, where n devices are grid-following 
converters or conventional energy storage converters; m 
devices are grid-forming devices. In the small-signal sta- 
bility analysis, grid-forming devices can be equivalent to 
an ideal voltage source with the series inductor, whose 
equivalent admittance is B1BR (j=1,..., m). Moreover, con- 
sider decoupling the actual operating conditions of the re- 
maining n devices from their dynamic. Thus, the inverse 
matrix of the system closed-loop transfer function is 

M =Y pr (s) +Y, a (s) 


where Y,,, (s) = diag (U? |P) Boa @¢(s) 
了 me (s)=diag (Zr (s)) 


where Bea represents the Thevenin equivalent admittance 
matrix only containing the n buses and considering the 
equivalent admittance Brsr;; Yrsri(s) is the admittance 
model of the i-th converter, i=1,...,1; P; and U; are power 
output and terminal voltage of the i-th converter. 


(1) 


An equivalent system matrix M; can be constructed 
in (2), whose eigen-subspace is approximate to its of the 
actual system matrix M since M can be treated as a 
perturbation of M; based on the perturbed theory of eigen- 


subspace. Moreover, the eigen-subspace R(X, and 


R(Y,) of the 
diag(U?/P;) Brea ® h with respect to its own eigenvalues 


equivalent admittance matrix 
A; is the same as matrix Mx. Therefore, using R(X,) and 
R(Y,) (KE {1,...,n}), a first-order approximation of the 
generalized eigen-matrix can be obtained as (3). Let E; 

be the closed-loop transfer function of the eigen-subsys- 
tem, the eigen-subsystem can approximate the stability of 
the actual system in the corresponding mode. The physical 
meaning of eigen-subsystems is a decoupled subsystem of 
the actual system since it physically retains the dynamic 


characteristics of converter and networks 
M, =1, 8G, (s)+Y,., (s) 


L, =Y} MX, =G, (s)+4 -¢(s) 
G, (s)=Y” | diag ((U}/P)Y (s))] X, 


=$ pu (U? (PY (s), Pri = Xa Yri 


where X= XIQ7， Vi= yh, the column vectors of X; 
and Y, are the right and left eigenvectors of 
diag(U?/Pi)Brea ® D with respect to As; xri and yri are the 
i-th elements of x; and yx, respectively. 


(2) 
(3) 


where 


On this basis, in small-signal stability mode, the sys- 
tem strength depends on the gSCR of the power grid, 
which is defined as the equivalent admittance correspond- 
ing to the eigen-subspace under the 1-th mode of interest: 


gSCR=y;" diag (U? /P)B „x 
=A (diag(U? /P)B a) 


represents solving for the minimum posi- 


(4) 


where Afin 
tive eigenvalue of matrices. 

When the dynamic G, (s) of the equivalent con- 
verter remains unchanged, the smaller the gSCR, the 
weaker the system strength. Based on this monotonic 
characteristic, the minimum gSCR required for the equiv- 
alent converter to reach a critical stable state can be deter- 
mined as (5). SCRo is precisely equal to the critical value 
of the gSCR(CgSCRo). 


SCR, = arg {det (G, (s,)+gSCR-E(s, )) = 0} (5) 


gSCR 
According to (6) and (7), it can be determined 
whether the system is stable and whether the system meets 
the stability margin requirements. 


gSCR> CeSCR, = SCR, (6) 
SCR_SC 
pes SO 100% 2 B% (7) 


The effectiveness of the system strength quantifica- 
tion method is verified in cases. 


